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ABSTRACT 


A  Fortran  program  is  presented  which  computes  a  normalized  equivalent 
source  strength  for  the  Class  I  and  approximates  the  Class  II  Flextentional 
Underwater  Acoustic  Transducer  Shells.  The  program  utilization,  the  input 
data  and  the  output  data  formats  are  described  in  detail. 
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INTRODUCTION 


This  report  describes  a  Fortran  program  of  an  equivalent  source  strength 

1 

model  for  the  Class  I  and  II  (Reference  1)  Flextensional  Underwater  Acoustic 
Transducer  Shells.  Pictures  of  these  two  shell  configurations  are  shown  in 
Figures  1  and  2. 

An  approximate  analytical  model  for  the  Class  I  type  of  shell  design 
has  been  developed  (References  2  and  3.  The  math  model  described  in  Refer¬ 
ences  2  and  3  utiiizes  a  numerical  model  in  order  to  obtain  the  near-  and 
far-field  acoustic  pressures  for  a  specified  displacement  on  the  shell  and 
associated  frequency  of  radiation.  This  model  is  sufficient  for  engineering 
applications  up  through  the  fundamental  mode  of  vibration  provided  that  a 
so-called  forbidden  frequency  is  not  encountered  (Reference  4). 

It  is  not  uncommon  for  the  effective  ka  of  a  flextensional  shell  to 
be  less  than  one  thus  yielding  a  relatively  omni-directional  far-field 
radiation  pattern.  If  for  such  design  situations  an  equivalent  source 
strength  model  was  available,  an  approximation  to  the  far-field  pressures 
at  a  defined  frequency  and  shell  displacement  would  be  possible.  In  addition, 
the  equivalent  source  strength  model  would  be  more  economical  to  utilize, 
computer  cost,  etc.,  and  would  yield  important  preliminary  design  information. 
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FIGURE  1 


FIGURE  2 


PROGRAM  UTILIZATION 


For  a  given  transducer  having  the  general  cent  iguration  depicted  in  F*-B“ 
ure  3,  all  geometric  quantities  pertinent  to  the  solution  of  the  free 
vibration  eigenvalue  problem  and  1 1  to  computation  of  Lite  source  strength 
may  be  determined  by  specifying  any  <<no  of  several  dimensionless  parameters 
and  the  opening  angle,  ^ .  The  parameters  L/D  and  0^  are  the  independent 
variables  in  the  accompanying  FORTRAN  program;  L/D  being  selected  since  it 
is  also  relevant  with  respect  to  the  size  of  the  piezoelectric  stack  re¬ 
quired  to  drive  the  transducer. 

The  solution  of  the  approximate  dynamic  model  involves  computing  the 
first  eigenvalue  and  eigenvector  of  a  2*  NP-2  square  matrix,  where  NP  is 
the  number  of  control  points  input  for  generation  of  the  discrete  system 
equations  of  motion.  It  is  suggested  that,  for  best  numerical  results, 

NP  should  lie  in  the  range  10  _  NP  ^  15. 

The  units  used  are  immaterial  since  all  input  -  output  quantities 
are  dimensionless  ratios  with  the  exception  of  the  opening  angle  which  is 
in  degrees. 
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INPUT  DATA 


On  the  first  clara  card  enter  NN,  the  number  of  cases  to  be  run  con¬ 
secutively.  NN  is  input  in  FORMAT  (12). 

On  each  following  card  enter  the  quantities  1VC.NP,  IAL,  i'AU,  INC,  RK , 
and  RATIO  in  FORMAT  (U,  12,  dll,  2F4.1). 

IVC  controls  the  printout  of  the  components  of  the  eigenvector  normal 
ro  the  surface  of  the  transducer  and  the  normalized  radiating  areas.  These 
vectocs  are  printed  out  only  if  IVC  is  entered  as  a  non  zero  positive 
quantity . 

NP  is  the  number  of  control  points  to  be  defined  on  one  half  of  a 
typical  transducer  element  for  the  discrete  system  dynamic  model. 

IAL,  IAU ,  and  INC  correspond  to  the  first  opening  angle,  the  last 
opening  angle,  and  the  opening  angle  increment,  respectively,  desired  in 
the  output  for  a  given  RATIO.  These  quantities  have  units  of  degrees. 

RK  is  the  ratio  of  the  radius  oi  curvature  of  a  typical  transducer 
element  to  the  radius  of  gyration  of  a  typical  element  section  as  measured 
at  the  guided  end,  The  dynamic  model  is  insensitive  to  this  parameter  for 
RK  ->  20.  RK  is  normally  input  as  any  large  positive  number  for  thin  transducer 
elements.  RATIO  -  L/D. 

PRINT-OUT  DATA 

The  input  quantities  RATIO  and  NP  appear  as  header  information  for  each 
case.  The  output  quantities  THETA,  RADIUS/LENGTH,  HEiGHT/EENGTH,  DiSPL. 

RATIO,  WIDTH  RATIO,  WIDTH  FACTOR,  COEFFICIENT  AND  NORM.  SOURCE  SIR-  ate 
printed  out  for  each  opening  angle  where: 

THETA  =  0 

o 

RADIUS/LENGTH  «  R/L 

HEIGHT/LENGTH  -  H/L 
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DISPI..  KA1 10  is  die  ratio  of  the  displacement  measured  at  the  center 
of  the  circular  arc  formed  by  a  typical  transducer  element  to  the  corresponding 
horizontal  displacement  of  the  guided  end. 

WIDTH  RATIO  is  the  ratio  of  the  width  of  a  typical  element  as  measured 
at  the  center  of  of  circular  arc  to  the  width  of  the  element  as  measured  at 
the  guided  end. 

WIDTH  FACTOR  -  element  width  at  center  of  arc/ (1/2  arc  length  of  element 
*  sin  (Pl/N)),  where  N  is  the  number  of  transducer  elements. 

COEFFICIENT  is  the  frequency  coefficient  i.e.  that  number  which,  when 
multiplied  by  SQRT  (E*I/m) / (L*L) ,  yields  the  natural  frequency  in  rad. /sec. 
where  E  is  the  modulus  of  elasticity,  I  and  m  are  the  area  moment  of  inertia 
and  mass  per  unit  length  of  a  typical  elemenc  as  measured  at  the  guided  end, 
and  L  is  as  shown  In  Figure  3. 

NORM.  SOURCE  STR  is  the  source  strength  of  the  transducer  normalized 
with  respect  to  ( (PI*D) A*2) *F*  DELTA  where  F  is  the  frequency  of  excitation 
in  cycles/sec-,  DELTA  is  the  displacement  amplitude  of  the  guided  end,  and  D  is 
as  in  Figure  ,3.  Negative  values  of  this  quantity  indicate  that  the  net 
volume  displacement  is  negative  when  the  ends  of  the  transducer  are  displaced 
in  such  a  manner  that  the  transducer  is  elongated. 

If  IVC  •-  0,  the  vectors  NORMAL  VECTOR  and  NORMALIZED  AREAS  are  also 
printed  out-  NORMAL  VECTOR  is  the  vector  of  eigenvector  components  normal 
to  the  transducer  surface  normalized  with  respect  to  Li.e  eigenvector  component 
at  the  guided  end,  NORMALIZED  AREAS  is  the  vector  of  radiating  areas  associated 
with  the  discrete  system  control  points  normalized  with  respect  to  the  area 
of  the  circular  end  section. 
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FINAL  REMARKS 


Ihe  model  presented  m  this  report  is  an  approximate  model  ii net- 
several  simplifying  assumpt  ions  were  necessary.  A  list  of  the:  most 
important  assumpt Ions  made  and  where  deemed  appropriate  the  aftf'icipated 
results  ol  having  to  make  such  assumptions  axe  os  follows: 

a)  1c  is  assumed  that  the  external  geometry  of  the  ■shells  can 

be  approximated  by  cylindrical  end  sections  (actual  lengch  is  ‘riot  considered 
since  It  is  assumed  dial  k«i<  <  l)  and  circular  curved  beams  exhibiting  a 
varying  taper  characteristic. 

b)  The  cylindrical  end  sections  do  hot  radiate  except  on  the 

•  ■'  ■'  ,  ■  ■  y.':  '  ■■  ..  ;•  ’ 

ends.  If  the  radial  thickness  of  the  cylindrical  end  sections  is  sniail 
compared  to  the  thickness  of  the  curved  beams  then  most  likely  these  end 
sections  will  contribute  significant! y  to  the  overall  radiation  field. 

c)  Raidation  through  the  slits  between  the  curved  beams  is 
jassumed  to  be  jinsigni  t  leant  . 

|  d)  It  is  assumed  that  the  boundary  condition  l:oi  die  curved 

ibe3ms  is  clamped -gui Jed .  Most  Class  I  shells  will  approximate  this  bound¬ 
ary  condition  as  will  relatively  thick  Class  II  shell  designs. 

e)  The  present  program  dues  not  allow  the  option  of  specifying 
the  actual  curved  beam  end  mass  ff  the  end  masses  become  large  as  com¬ 
pared  to  the  total  mass  of  the  curved  beams,  as  is  quice  possible  f.ot  the 
Class  II  shell,  then  the  eigenvector  associated  with  the  shells  fundamental 
frequency  of  vibration  approaches  the  quasi-static  solution.  This  would 
result  in  a  reduction  in  die  predicted  far-field  pressures  of  approximately 
25%. 

The  program  desorbed  hoe  in  has  been  checked  by  comparing  the  predicted 
and  measured  far -f ie Id  pressures  tu:  a  transducer  of  known  performance ,  The 
predicted  pressures  were  approximately  M.5  db  above  die  measured  data  In 
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1  iglit  of  the  intended  use  for  the  model  described  herein  the  above  variation 

is  acceptable.  In  addition,  since  it  is  entirely  possible  to  design  a 

I  1  extens lonai  transducer  (.for  ka  .' )  with  0  0,  this  program  should  be 

s 

utilized  early  m  the  design  process  in  order  to  elimate  this  type  of 
design  error. 
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APPENDIX  A 


Listing  of  the  Foctran  Program 


C  THIS  PROGRAM  COMPUTES  THE  SOURCE  STRENGTH,  FREQUENCY  COEFFICIENT, 

C  AND  A  NUMBER  OF  GEOMETRIC  PARAMETERS  FOR  A  CLASS  OF  FLEXTENSIONAL 

C  ACOUSTIC  TRANSDUCERS  COMPOSED  OF  AN  ARRAY  OF  TAPERED  CURVED  BEAMS 
C  (CIRCULAR  ARCS)  FORMING  A  SURFACE  OK  REVOLUTION  AND  EXHIBITING 
C  SHORT  CYLINDRICAL  END  SECTIONS 

C  USAGE-  THE  FOLLOWING  PARAMETERS  MUST  APPEAR  ON  THE  DATA  CARD 
C  FOR  EACH  CASE-  I  VC,  NP,  IAL,  I.AU,  INC,  RK,  RATIO 

C  DESCRIPTION  OF  INPUT  PARAMETERS 

C  NN= NUMBERS  OF  DATA  SETS 

C  TVC=I  IF  A  PRINTOUT  OK  THE  NORMALIZED  AREAS  AND  NORMALIZED  NORMAL 

C  VECTORS  IS  DESIRED (THE  TANGENTIAL  COMPONENTS  OF  THE  EIGENVECTOR  ARE 

C  DESTROYED  IN  THE  NORMALIZATION) ,  IVC=0  OTHERWISE 
C  NP=THE  NUMBER  OF  CONTROL  POINTS  IN  HALF  OF  THE  BEAM 

C  TAT -FIRST  OPENING  ANGLE  DESIRED  IN  THE  OUTPUT 

C  TAU=LAST  OPENING  ANGLE  DESIRED  IN  THE  OUTPUT 
C  TNC=OPF,NING  ANGLE  INCREMENT 

C  RX=BEAM  RADI US/ RADIUS  OF  GYRATION  OF  A  TYPCIAL  ELEMENT  EVALUATED 
C  AT  THE  GUIDED  END 

C  RATIO*CHORD  SUBTENDED  BY  THE  OPENING  ANGLE/DIAMETER  OF  CYLINDRICAL 
C  END  SECTION  AT  THE  GUIDED  END 

C  DESCRIPTION  OF  OUTPUT  PARAMETERS 

C  T H ET A = 0 P E N T N G  ANGLE  SUBTENDED  BY  A  TYPICAL  ELEMENT  OF  THE  ARRAY 

C  RADIUS/LENGTH®3EAM  RADIUS/CKORD  LENGTH 

C  HEIGHT/LENGTH*TRANSDUCER  HEIGHT/ CHORD  LENGTH 

C  DISPL.  RATIO® RATIO  OF  NORMAL  DISPLACEMENT  AT  THE  BEAM  CENTER  TO  THE 
C  DISPLACEMENT  OF  THE  GUTDED  END 

C  WIDTH  RATT0=RA7T0  OF  BEAM  WIDTH  AT  THE  CENTER  OF  THE  BEAM  TO  THE 
C  BEAM  WIDTH  AT  THE  GUIDED  END 

C  WIDTH  FAC70R=MAXTMUM  BEAM  WIDTH/ (ARC  LENGTH  OF  HALF  BEAM*S  IN(PI/ 

C  1N)),N=NUMBER  OF  BEAMS  IN  THE  ARRAY 
C  COEFFICIENT® FREQUENCY  COEFFICIENT 

C  NORM.  SOURCE  STR.®SOURCE  STRENGTH/ ( ( (PI*D) **2) *F*DEL  TA) ,  D®END 

C  DIAMETER,  F= FREQUENCY  IN  CY./SEC. ,  DELTA®DI SPLACEMENT  OF  GUIDED  END 

C  NORMAL  VECTOR-NORMAL  EIGENVECTOR  COMPONENTS  NORMALIZED  ON  THE 
C  DISPLACEMENT  OF  THE  GUIDED  END 

C  NORMALIZED  AREAS=RADIATING  AREAS  ASSOCIATED  WITH  EACH  DISCRETE  CONTROL 

C  POINT  NORMALIZED  WITH  RESPECT  TO  THE  AREA  OF  THE  CYLINDRICAL  END 
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MATN’  (c.'ti^inucd') 


DIMENSION’  A('^.K)  %RM(uO>  .KI (GO) ,RN(20>  ,RS(30) 

COMMON  N’.M,  MODI'.,  A.  VU 
PI  =  3.  1.;‘,:AVW 
M=  I 

MODE®] 

READ(1  .  i02):,'N 
102  FOICwVf'I/} 

DC  ] 


READ  ( t  ,  1 03  )  I  VC ,  M? ,  I  AT. ,  I AU ,  INC .  RK ,  RATI  0 

103  FORMAT (Ti » 12 , J 13 ,2F4 . 1) 

N*2*NP-2 

WIT  E  (  3 . 1 00  >  RATI  0 ,  NP 

100  FORMAT ('T' //////& 7X,' LENGTH/ END  DIA.-'F7.4/47X,  'NO.  OF  CONTROL  PTS 
1.  IN  HALF  BEAM® ' T2////8X , 'THETA' ,  3X , ' RADIUS/ LENGTH* ,3X, 'HEIGHT/LEN 
2GTH'  ,3X.  'DISPL.  RATIO5  ,3a,  'WIDTH  RATIO'  ,3X, 'WIDTH  FACTOR'  ,3X,  ’COEF 
3FJCIENT' ,3X,"N0RM.  SOURCE  SIR. ') 

DO  1  LA-IAL , TAU, INC 
0A-LA*?T'180. 

RR-SIN(OA/2.) /RATIO 

CALL  COEF F(N , N*P , OA , RR ,  SI , PHI  ,RM,RI, RN , RS) 

CALL  STIFF ( N , NP , RK , SI , PHI , A , RM , RI , RN) 

CALL  GIVENS 

VU(  1 )  *  SORT (VI! ( !) > *COS ( PHI / 2 .  ) *<S  IN(0 A/ 2 .  ) **2 ) / < STN< PHI  / 2 . .) **2) 

DO  2  I-1,K 

A(T,1)-A(t,l)/SQRT(RM(I» 

2  CONTINUE 

qsN“A(1  ,  i  Va?s(a(  1 .  i )} 

DO  3  1*2. N, 2 
K-l/2+1 

QSN-QSI+t-AC  .  ln.  *RS(X)  /ABS(A(1 , 1)) 

3  CONTINUE 
TF(A(1 , 1/ .4,5 

4  QSN*-f<SN 

5  RL*1. / (2, *SIN(0A/2.  )) 

THL- (RR+ 1 . -COS (OA/ 2 .  > ) / SIN (OA/ 2 .  ) 

DR-A (N , 1 ) /A (1,1) 

WR*RI(N?) 

WF=4.*(1. -COS(OA/2. ) >STN(OA/2.  ) /RATIO) /OA 
OA=OA*1RO. /PI 

WRITE 0 . 1 0 1 )  OA .  P.L ,7HL , DR ,  WR , WV ,  VO(  1 )  , QSN 
101  FORMAT (/ /F13. 1 , 2E16.  6 , El 5. 5 , E14 . 3 , E15. 5 , E14. 5 , El S.  6) 

TF(IVC) 1 , 1 .6 

6  IF(A(l,l)?7f7f8 

7  A(l,l)--A(lfl) 

DO  9-T-2.N.2 
A(I,1)— A(T.l) 

(15; 


MAI1.  (continued) 


9  CONTINUE 
8  DO  10  T=2 ,N,2 
K-T/2+1 

A(K,1)-A(T,1)'A(1,1) 

10  CONTINUE 
A(1 ,1)=1 . 

WRITE (3 , 10A; 

10A  FORKAT(//34X, 'NORMAL  VECTOR' ,34X, 'NORMALIZED  AREAS'//) 
WRTTE(3, J05)(A(I,1) ,RS(I) ,I=1,N?) 

105  FORMA?(35X,E12. 5,36X,E12. 5) 

1  CONTINUE 
END 
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COEl-'F 


SUBROUTINE  COEFF(N ,NP ,OA ,RR , SI ,PHI ,RM ,RI ,RN ,RS) 

DIMENSION  RM(60) ,RI(30) ,RN(30) ,RS(30) 

Pl=>3. 1 415927 
SI* (PI-OA) /2 . 

PHI*OA/ (2 . * (NP-1.  ) ) 

OA2=OA/2. 

PHI2*PHl/2. 

D«COS(OA/2.)-RR 

RI(1)-1. 

RN  (].)*"  1 .  RS(1)«1 

DO  1  J=2 ,NP 

RI (J)= (COS (OA2- (J-l. )*PHI)-D)/RR 
RN(J)=2.*RI(J)*R1(J-1)/(RI(J)+R1(J-1)) 

RS  (J )“2.  *  ( 2 .  *SIN  (P1II2 ) *COS (OA2- (J-l.  ) *PHI)-D*PHI) /  (RR*RR) 

1  CONTINUE 

RS(2)*2. * (2 . *SIN(3 .  *PHI/4.  )*COS (OA/2.  -3.  *PHI/ )-D*3.*PHI/ 2. )/  (RR* 
1RR) 

RS (NP)*2. *(SIN(PHI2)-D*PHI2)/(RR*RR) 

RM ( 1 ) * 1 . 

RM(N)=RI(NP)/2. 

M-N-2 

DO  2  J*=2  ,M,  2 
RM(J)*RI (J/2+1) 

RM(J+1)*RM(J) 

2  CONTINUE 
RETURN 
END 


an 


STIFF 


SUBROUTINE  Sl  l  FF(N,Nl*,KK,SI , PU1  ,A,RM,R1,RN) 

DIMENSION  A(60,60> ,RM(60) ,RN(30) ,Rl(30) 

C*SIN(PHi/2 . ) / COS (PHI/ 2 . ) 

C2*SIN(SI)/COS(SI) 

P*  (  2  .  *RK  *  S 1 N  ( I’ll  I  /  2  . ) )  *  4  2 
DO  1  1=1, N 
DO  1  .1=1,  N 
A ( 1  ,  l)=0.0 
t  CONTINUE 

A  ( 1. ,  !  )  =  (  (C*  {C  -‘UN  (2)  *I'~  (2 .  fRIU))*(C*C2-l.))*COS(SI)+((C+C2*RN 

1  ( 2)  *P+-(2  .  +  R r  (2) 1 *C* (C*C2-1 . ) )*SIN(SI) )  *COS (SI)  :  ‘t: 

AO  ,2WC*(C+C2)*RN(2)*P+2  .*(1  .+RI  (2)  )  4  (C*C2-1 . )  )*COS  (SI) 

A(l,3)**(  (C+C2) *RN  (  2)  *P-2  ,*C*  (C*C2-1 .  )  )  *COS  (SI) 

A(1 ,  A )  =-RI.  (  2) *  (C*C2-1 .  ) *COS  (SI) 

A(lt5>»RI(2>*C4(C*C2-l.)*COS(Sl) 

A (2 , 2) - (RN ( 2 )  +  RN ( 3) ) *P*C*C+2 . +4  . *RI (2)+RI (3) 

A 12 , 3)=C* ( (RN(2)-RN (3) ) *P-2 .+RI (3) ) 

A(2,4)«RfJO)*P*C.*C-2.*(RT  (2)+RI(3)) 

A(2,5W;*(RN(3)*P*2.*RI (2)) 

A(2,6)»Rl(3) 

A(2-,/)— RI(i)*C 

A(3,3)*(RN(2)+RN(3))*P+(2.+RI(3))*C*C 

A(3,A)—C*(RNO)*P*-2.*RI(3)) 

A(3 , s)--RN(3)*P 

A(3 , 6)=RI (3)*C 

A(3,7)— RL(3)*C*C 

KI-N-A 

K2=N-3 

K3-N-2 

K4=N-1 

Jl=NP-3 

J2-NP-2 

I3*NP-1 

A0U,K1)«(RN(J2)+KN(J3))*P*C*C+RI(J1)+4.*R1  (J2)+RI(J3) 

A(K1,K2)=C*  ( (ICN  ( J2)-RN  (J3) )  *P-R1  (Jl)+Rl(J3)  ) 
A(Kl,K3)!-RN(J3)*r*C*C-2.*(RI(J2)+RI(J3)) 

A(K1-,K4)«C*(RN(J3)*P+2.*RI(J2)) 

A(K1,N)«RI(J3) 

A(K2,K2)  =  (RN<  J2)+RN(J3) )*P+(Kl (Jl)+RI (J3) )*C*C 
A(K2,K3)—  C*(RN(J3)*P+2.*RI(J3>) 

A(K2,K4)— KN(.I3>*P 
A(K2,N)**Rl  ( J  3)  *C 

A(K3,K3)=* (RNl  J3 )  fRN(NP) )*P*C*C-fRI (J2)+4 . *R1  (J3)+2  ,*RI (NP) 

A(K3,K4)=C*  ( (RN  ( J3)-RN(NP)  )*I’-RI  (J2J4-2 .  *RI  (NP)  ) 
A(K3,N)“RN(NP)*P*C*C-2.*RI (J3)+RI (NP) ) 
A(K4tK4)-(RN(J3)+RN(NP))*P+(RI(J2)+2. *RI(NP))*C*C 

A(N',N)-RN(NP)*P4C*C-FRl(J3)+2-.*RI.(NprVfA(K4,N)— C*(RN(NP)*P+2.:*RI(NP)) 
JM-N-6  L 

DO  2  1*4, JM, 2 
DO  2  J  *  A , JM , 2 

i  F  ( I  -  J  )  2 , 3 , 2 
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S II  i'T  (cent  luued) 


3  K-I/2+1 

A( I , J)= (RN(K) 1RN  (K+l)  )  *p*c*C4-RI  (K-l)+4  .  *RI  (K)+RI  (K+l) 
A(I , J+1)=C*  ( (RN(K)-RN (K+l)  )  *P-RI (K-1)+RI  (K+l)  ) 

A(I , J  +2)  *RN  (K+l ) *P*C*t,-2 . *  (R1  (K) +RI  (K+l)  ) 
A(I,J+3)*C*(RN(K+1)*P+2.*RI(K)) 

A(I , J+4)=RI (K+l) 

A(I,J+5)=-RI(K+l)*C 

M*  1+1 

L=J+1 

A (M ,L) * (RN (K) +RN <K+1 ) ) AP+ ( RI (K-3 )  +RI (K+l) ) *C*C 
A (M ,  L+l )»*-C*  (RN  ( K+l  )  *P+2  .  * P I  (K+l )  ) 

A(M,L+2)*-RN(K+l) *P 
A(M,L+3)-RI(K+l)*C 
A (M , L+4 ) *-RI ( K+l ) *C  *  C 
2  CONTINUE  1 
1  DO  4  1*1, N 
DO  4  J*1,N 
A ( J , 1 ) *A ( 1 , J ) 

4  CONTINUE 
DO  5  1*1, N 
DO  5  J  * 1 , N 

A(I,J)*A(I,J)  /  SQR 1'  (  RM  ( I )  *  KM  ( .1 ) ) 

5  CONTINUE 
RETURN 
END 
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CIVENS 


SU BKOU'I  I NE  G  !  VENS 

DIMENSION  A(60,60) ,E(60,60) ,DIAG(60) ,VU(60) .VLCeoTTiDCeO) , 

,S(60)  ,C(60)  ,D(60)  ,IND(60)  ,U(60) 

COMMON  N  ,H .  MODE ,  A ,  VU 

EQC  t  V a  CNCIi  C P ,  PRODS)  ,  (TAU  ,  BETA)  ,  (T , SMALLD)  ,  (II , MATCH)  ,  (VL(1)  ,D(  1) 

l).{')U),S(l))  ,|A 

C  CAI.Oll.AT I:  NORM  OF  MATRIX 

ANORM  0. 

DO  1  l«|,N 

sn( i ) ~o . 

DO  l  J  » 1  ,  N 

l  ANORM- ANUKTH-A(  I  ,J)*A(l  ,J) 

A  1  KM-  "QIC-.  <  ANORM) 

C  GENERATE  IDENTITY  MATRIX 

ii-uoioooa  1,1000 

1000  DO  7  I»L,N 
DO  7  J - I , N 
IF( I-J )6, 1001 ,6 
100]  60 

GO  TO  7 

6  B(1  jj)«0. 

7  CONTINUE 

I 

C  PERFORM  ROTATIONS  TO  REDUCE  MATRIX  TO  JACOBI  FORM 

11  NN*N-2 

LF(MN)  100,21. ,! 003 

1003  DO  15  1=1 , NN 
11*1+2 

DO  1 5  J= II  . N 
T1«A(I t LEI) 

'  T2«AU,J) 

1  FT  T  2) 1 0OP ,15, 1004 

1004  CALL  ROTATE (T1 ,T2, SN ,CSN) 

DO  12  K*l ,N 

T2~-CSN*A(K ,  I  f  l)+SN*A(K,J) 

A(K,J)*CSN*A(K,J)-SN*A(K,1+1) 

12  A(K. l  +  l » —  r 2 
DO  ! 3  K-f ,N 

T2«CSN*A(I+I ,K)+SN*A(J,K) 

A0.K)*CSN#A(J,K)-SN*A(I+1,K) 
n  A(  f  l  ,K)*T2 

I  F (ft)  1 00 5,15,1005 

1005  DO  !4  K* I ,W 

T2--CS!v  B(K , !  •*  I  )+SN*B(K.J) 
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GIVENS  (continued) 


3(K,J>*CSN*B(K,J)-SN*B(K,I+1) 

14  B(K,l+l)=T2 

15  CONTINUE 

C  MOVE  JACOBI  ELEMENTS  AND  INITIALIZE  EIGENVALUE  BOUNDS 

21  DO  22  I- 1  ,N 
DIAG(1)-A(T,r) 

VU(T)-ANCRM 

22  VL(l)— ANORM 
DO  23  1=2, N 
SD(I*l)=Aa-l,l) 

23  Q(I-1)«SD(I-1)*SD(1-1) 

C  DETERMINE  SfCNS  OF  PRINCIPAL  MINORS 


TAU=0 

1=1 

35  MATCH-0 
T2=0. 

Tl-1. 

DO  2b  J= l.N 
P=D I AC ( J ) ■ TAU 
IF(T2) 1007,27, 1007 

1007  IF(T1) 1008, 28 , 1008 

1008  T-P*Tl-Q(J-l)*T2 
GO  TO  29 

27  IF(T1) 1009 ,30,30 

1009  Tl— 1. 

T— P 

GO  TO  29 

30  Tl-l. 

T=r 

GO  TO  29 

28  IF(Q(J-1))1010,30,1010 

1010  IF(T2)31 ,1011,1011 

1011  T— I, 

GO  TO  29 

31  T-l. 

C  COUNT1  AGREEMENTS  IN  SIGN 

29  IF(T1)32, 1012,1012 

1012  IF (T)34 ,33,33 

32  IF(T)33 ,34,34 

33  MATCH-MATCH-)- 1 

34  T2-T1/ABS(TI) 

26  T1-T/ABSCT1) 
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GIVENS  (continued) 


:  ESTABLISH  tighter  bounds  on  eigenvalues 

do  <»  i  K“I  ,  MODE 
I F (K*  N+MATCH) 1014,1014 ,42 

1014  IFtVUOO-TAOM  ,41,1015 

1015  VU(K)»TAU 
GO  TO  ..l 

42  TFtVUK)  -1AL)  1018,41 .41 

1018  V 1. 1 K  )  “TAl1 
M  CONTINUE 

45  I  F(VU(1 )  *  VL(  1 )  • 5. E-6)43 ,43 , 1016 

1016  iF(VUvO)  1017,44,1017 

1017  17 (ABS\  Y'L(  I ) /VU( I)-l,  V5.  E-6)43 ,43 ,44 

43  1  =  1+1 

TF(I-MODE')  1019,1019,51 

1019  GO  TO  ^5 

44  TAU-vVL'.  I)+VU(1));  . 

GO  10  35 

51  IF(M)'  *00,100, 1020 

1020  1 F \M -  MODE) 1021,1021,1022 

1022  M-M0I1E 

1021  DO  52  l«l,N 
DO  52  J= 1 , N 

52  A(I,J)=0. 

DO  53  1=1 ,M 
IF(I-l) 1029,54 ,1029 

1029  IF rvU(l)  •  VlKT-D-5. E-6) 55,55,1023 

1023  I F(VUC l ) ) 1024 , 54 , 1024 

1024  IF(ABS(VU(l- l)/VU(I)-l. )-5. E-6)55 ,55 ,54 

54  CSN-l. 

SN=0 

DO  56  J  =  I. , N 
IF(J* 1)1025,57,1025 

1025  CALL  ROTATE (Ti ,T2 , SN , CSN) 

S(J-  l)=SN 
C(J- l )=CSN 
D(  J-l)=7  !*CSN+T2*SN 

57  T 1  ■  *(D  I  AG (  J)  *  VU ( 1  )  )  *CSN-LLiA"SN 
T2*SD(J) 

56  BETA=SD(J) *CSN 
D(N)"T 1 
DO  58  J*'l  ,N 

58  IKD'(J)»0 

55  SMALLD«ANORM 
DO  60  J  ~  i  ,  N 

t!?(lNI)(  J ) *  I  )  1026,60,60 

1026  I F( AP.S ( SMA1.LD) -  ABS (D(J) ) ) 60 , 60 , 1027 

1027  SMALL  0*4)  (J) 


/no\ 


GIVENS  (continued) 


NN~J 

60  CONTINUE 
INIHNN';-] 

IF(NN- I } 1028,53 , 1028 
1028  DO  61  K- 2 , NN 
! 1=XN+1 *K 

A(  11+1,1)  ~C(  1 1  )  •’ PRODS 

61  PRODS--  PROPS ’*S(  11) 

53  Av l , 1 ) -PRODS 

C  FORM  MAI U IX  PRODUCT  OF  ROTATION  MATRIX  WITH  JACOBI  VECTOR  MATRIX 

DO  6b  .1-1  ,M 
DO  6  7  K~  1  , N 
67  U ( K) =A ( K , J ) 

DO  6b  T- 1  , N 
A(I  , J)-0. 

DO  06  K - 1 , N 

66  Ai  1  ,J)"Bv  I  ,K)  *IKK)+A(I  ,J) 

100  RETURN 
END 


S U BRO IJTINi:  R01ATE(A.B,C,D) 
E=SQRT  vA,:A+B*B) 

D»A/E 

C-B/E 

RETURN 

END 
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APPENDIX  B 

Listing  of  a  Typical  Output  Data  Set 
TYPICAL  OUTPUT 
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TYPICAL  OUTPUT  (continued) 
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